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1 Introduction 


Supersymmetry seems to be a necessary ingredient of a quantum theory of gravity. Many of the 
possible extensions of the Standard Model beyond the presently known energy range are based 
on supersymmetry, ft is generally assumed that the scale where supersymmetry becomes man¬ 
ifest is near to the presently explored electroweak scale and that the supersymmetry breaking 
is spontaneous. An attractive possibility for spontaneous supersymmetry breaking is to exploit 
non-perturbative mechanisms in supersymmetric gauge theories. This is the basis of a strong 
theoretical interest for investigating supersymmetry non-perturbatively. 

The motivation to investigate non-perturbative features of supersymmetric gauge theories 
is partly coming from the desire to understand relativistic quantum held theories better in 
general: the supersymmetric points in the parameter space of all quantum held theories are 
very special since they correspond to situations of a high degree of symmetry. The basic work 
of Seiberg and Witten JT|] and other related papers showed that there is a possibility to ap¬ 
proach non-perturbative questions in four dimensional quantum held theories by starting from 
exact solutions in some highly symmetric points and treat the symmetry breaking as a small 
perturbation. Beyond this, the knowledge of non-perturbative dynamics in supersymmetric 
quantum held theories can also be helpful in understanding the greatest puzzle of the standard 
model, with or without supersymmetric extensions, namely the existence of a large number of 
seemingly free parameters. As we know from QCD, strong interactions in non-abelian gauge 
theories are capable to reproduce from a small number of input parameters a large number of 
dynamically generated parameters for quantities characterizing bound states. This is a possi¬ 
ble solution also for the parameters of the standard model if new strong interactions are active 
beyond the electroweak symmetry breaking scale. 

The simplest supersymmetric gauge theory is the supersymmetric extension of Yang-Mills 
theory. It is the gauge theory of a massless Majorana fermion, called “gaugino”, in the adjoint 
representation of the gauge group. The Euclidean action density of a gauge theory in the 
adjoint representation can be written as 

£ = J W, + i (®„A)“ + m 5 TA“ . (1) 

Here F^ u denotes the held strength tensor and A a is the Grassmannian fermion held, both with 
the adjoint representation index a. nig is the gaugino mass which has to be set equal to zero 
for supersymmetry. For a Majorana fermion A a and A a are not independent but satisfy 

A = X T C , (2) 

with C the charge conjugation Dirac matrix. This definition is based on the analytic continu¬ 
ation of Green’s functions from Minkowski to Euclidean space [|J. 

The held theory defined by the Euclidean action ([I]) has for nig = 0 a “supersymmetry” 
with respect to the infinitesimal transformations with a Grassmannian parameter e: 

5AI = 2ze 7/ X , 5A“ = -a^F^e . 
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(3) 


It is easy to show that the change of the action density is a total derivative: 


&£ = e = e d idn 

where j p = and the “supercurrent” S p is defined as 

(4) 

S, = -F^a pTlil X a . 

(5) 


Note that this is a Grassmannian Majorana current having both a four-vector and a Dirac 
spinor index (this latter is not explicitly shown here). 

The existence of a supersymmetry in the above simple gauge field theory is at first sight 
surprising. It can be better understood in the general framework of supersymmetric field 
theories based on “superfields” [3|. In the Wess-Zumino gauge the action of Yang-Mills theory 
with N — 1 supersymmetry is conventionally given as 

J d 4 xd 2 9Tr(W a W a ) 


= J d i xTr^F lu ,F'" + l -F^F^ -iXa p (D p X) + i{D li X)a /i X + D 2 } , (6) 

where the first line is written in terms of the spinorial field strength superfield W(x, 9 , 9) a which 
depends on the four Minkowski-space coordinates x p , /i — 0,1, 2, 3 and the anticommuting 
Weyl-spinor variables 9 a ,9a (a, a = 1,2). After performing the Grassmannian integration on 
9, one obtains the second form in terms of the component fields which are in this notation 
represented by Lie algebra elements. For instance, the field strength tensor F pv and its dual 
are defined as 

F^{x) = -igF“ v (x)T a , F pv = ^e p „ pa F pa . (7) 

A, A in (Up represent the Weyl components of the gaugino field. In the equation ([]])-(||) we 
denoted the corresponding Dirac spinors by the same letter A. In the present paper we shall 
use, with very few exceptions, Dirac spinors therefore this will not give rise to confusion. 

The action in (|p includes a 0-term, therefore it is natural to introduce the complex coupling 


0 47ri 

T= 2i Ar lfl 


( 8 ) 


and then, with arbitrary 0, the N — 1 SYM action becomes: 


r j d 4 xd 2 9Tr(W a W a )} 


d A x Tr 


-~F pv F pv - iXa^D^X) + i(D p X)d p X + ZY 


0 

167T 2 


d A x Tr 




(9) 


Performing the trivial Gaussian integration over the auxiliary field D } going to Euclidean space 
and setting the ©-parameter to zero one obtains the action in the form ( lj). 
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The supersymmetry transformations in (|3|) are relating the bosonic gauge field to the 
fermionic gaugino field. The above form of the transformation is “on-shell” because the aux¬ 
iliary field D is eliminated by the equation of motion D = 0 and it is “non-linear” due to the 
Wess-Zumino gauge fixing || . The realization of the supersymmetry in quantum field theory is 
not trivial because of the supersymmetry breaking introduced by the gauge fixing (for a recent 
discussion see |Q). Generally speaking one expects the existence of a renormalized supercurrent 
operator Sr^ satisfying the Ward-Takahashi-type identity 


duS R n = 2 m R XR , 


X = - 2 f, 


a 

tU'VfJ'V 


V , XR = Z x xj 


( 10 ) 


ffere rn R is the renormalized gaugino mass and Z x is an appropriate multiplicative renormal¬ 
ization factor. The consequences of supersymmetry at rriR = 0 can be obtained by considering 


different matrix elements of (|T0|) . The non-vanishing right hand side of the above equation 
describes the “soft breaking” of supersymmetry due to non-zero gaugino mass. 


1.1 Open questions of SYM dynamics 

On the basis of its similarity to QCD one can assume that the basic features of SYM dynamics 
are similar to QCD: confinement of the coloured degrees of freedom and spontaneous chiral 
symmetry breaking. As in QCD, a central feature of low-energy dynamics is the realization 
of the global chiral symmetry. There is only a single Majorana adjoint “flavour” therefore the 
global chiral symmetry of N — 1 SYM is U(1)a which coincides with the so called R-symmetry 
generated by the transformations 

0' a = e*e a , 0' & = e-*0 & . (11) 


This is equivalent to 


A' = e iip \ c 


A', = e~ i<p Xr 


A' = e~ i<pl5 X , A' = Ae"^ 75 


( 12 ) 


Here A a , denote Weyl components and A, A without spinor indices are the Dirac spinor 
fields. 

The U(l) A-symmetry is anomalous. For definiteness, let us consider in what follows SU(N C ) 
as gauge group when the corresponding axial current J M = Ay^qsA satisfies 


<9% 


N a 2 

C J {ivper rpr rpr 

32t r 2 ,MJ p ° ' 


(13) 


Here g denotes the gauge coupling. The anomaly leaves a Z 2 n c subgroup of U( 1)a unbroken. 
This can be seen, for instance, by noting that the transformations in (|I2|) are equivalent to 


rrig —> rrige 2 * ¥ ’ 75 , 0 —> 0 — 2 N c ip , 


(14) 
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where 0 is the ©-parameter of gauge dynamics. Since 0 is periodic with period 2ir, for nig = 0 
the {7(1 )a symmetry is unbroken if 

n 7T 

V = Vk = j^-, (k — 0,1,..., 2N C — 1) . (15) 

For this statement it is essential that the topological charge is integer. 

The discrete global chiral symmetry Z 2 n c is expected to be spontaneously broken by the 
non-zero gaugino condensate (XX) ^ 0 to defined by {<p 0 , Pn ( } (note that A —> —A corre¬ 
sponding to ip n c is a rotation). The consequence of this spontaneous chiral symmetry breaking 
pattern is the existence of a first order phase transition at zero gaugino mass rrig = 0. For 
instance, in case of N c = 2 there exist two degenerate ground states with opposite signs of the 
gaugino condensate. The symmetry breaking is linear in nig, therefore the two ground states 
are exchanged at nig = 0 and there is a first order phase transition. 

For larger number of colours (N c > 2) the phase structure is more involved. As an example, 
N c — 3 gives rise to three degenerate ground states for the first order phase transition at the 
supersymmetry point. (For a first numerical study of SYM with SU( 3) see [||.) 

There are analytical predictions for the magnitude of the gaugino condensate based on 
instanton calculations. On general grounds the result has to be proportional to A 3 where A is 
the dynamical scale developed by dimensional transmutation which is defined by 

A = ^-i/2/M/O 3 , & = ^ . (16) 

Here /j is the scale belonging to the coupling g(fi). As usual, /3 0 denotes the first coefficient of 
the /3-function. In terms of A 3 the magnitude of the gaugino condensate can be written as 

(A“A““) = CA 3 e 2nik/Nc . (17) 


The phase factor depending on the integer k refers to the different ground states defined in 
(]T5|) . The proportionality factor C depends, of course, on the renormalization scheme where 
A is defined. So called “weak coupling” instanton calculations I, 0-1 imply that we have for 
SU(N C ) gauge group C = 327T 2 in the dimensional reduction scheme A = A or [0- Another 
way of calculation using “strong coupling” gives, however, by a factor 2/((N c — 1)!(3Y C — 
l)) 1//JVc different result ]TD], [IT], [H|, |HSfl. (For a critique of the second method see ||14|| .) This 
discrepancy between the “weak coupling” and “strong coupling” results could perhaps be due 


to the existence of an additional chiral symmetric vacuum state [15 


The predicted magnitude of the gaugino condensate (EZP can be checked in lattice Monte 
Carlo simulations. For this it is convenient to switch to the lattice A-parameter A lat- First 


one can use 


a W a I 75 = exp{-l/18} 

and then for the lattice action which will be introduced in sections 0 and PI one has |[I6 

1 


( 18 ) 


a am/V.4t = exp 



N r n„ 

N C P + 
















P 3 = 0.0066960(1). 


N 

A) = ^( n - 2 n a ) , P = 0.0849780(1) , 


(19) 


Here n a is the number of Majorana fermions in the adjoint representation, that is for SYM we 
have to set n a — 1 . 

An important and interesting question is whether supersymmetry can be broken sponta¬ 
neously or not. In pure SYM theory, without additional matter supermultiplets, this is not 
expected to occur. An argument for this is given by the non-zero value of the Witten index Q 


IV = Tr (-1) F = n boson -n f 


ermion 5 


( 20 ) 


which is equal to the difference of the number of bosonic minus fermionic states with zero 
energy. It is supposed not to change with the parameters of the theory. For SYM theory we 
have Wsym — N c , therefore there is no spontaneous supersymmetry breaking. (The N c ground 
states discussed above correspond to this Witten index.) 

I 11 analogy with QCD, one expects that the spectrum of the SYM theory consists of colourless 
bound states formed out of the fundamental excitations, namely gluons and gluinos. (In this 
context we shall use the name “gluino” instead of the more general term “gaugino”.) In the 
supersymmetric point at zero gluino mass these bound states should be organized in degenerate 
supersymmetry multiplets. For the description of lowest energy bound states one can try to 
use an effective field theory in terms of suitably chosen colourless composite operators. 

For N — 1 SYM a low energy effective action was constructed by Veneziano and Yankielowicz 
(VY) [|1|]. The composite operator appearing in the VY effective action is a chiral supermul- 


tiplet S containing as component fields the expressions for the anomalies ]T9 


S = A(y) + V20My) + 0 2 F(y) , 


( 21 ) 


where y = x + idad. The scalar component of S is proportional to the gluino bilinear 


A oc X Q X C 


( 22 ) 


The other components contain gluino-gluino and gluino-gluon combinations. Therefore, as far 
as a constituent picture is applicable to the bound states formed by strong interactions, the 
particle content of the lowest supersymmetry multiplet is: a pseudoscalar gluino-gluino bound 
state, a Majorana spinor gluon-gluino bound state and a scalar gluino-gluino bound state. In 
terms of S the VY effective action has the form 

S VY = ^ J d 4 x d 2 9 d 2 9 (5+ S) 1/3 + 7 [J d 4 x d 2 9 (S log j - S) + h.c.] . (23) 

Here a and 7 are positive constants and A is the mass parameter for the asymptotically free 
coupling which also appears in (pT|). 

An interesting question is how the spectrum of glueballs, gluinoballs and gluino-glueballs 
is influenced by the soft supersymmetry breaking due to a non-zero gluino mass rrig A 0. For 
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small rrig it is possible to derive the coefficients of the terms linear in nrig in the mass formulas 
pOll: 


Mfj 

Mx 

M k 


N c aA + 
N c aA T 
N c aA + 


407T 2 

\ m ~g\ 

3A 

487T 2 

4 

\ m ~g\ 

3A 

567T 2 

4 

\ m g\ 


3 N c 


(24) 


Here fj , \ and / 0 denote the bound state with spin 0~, \ and 0 + , respectively. The range of 
applicability of the linear mass formulas in (]24|) is not known. 

The main assumption needed to derive the VY effective action (|23|) is the choice of the 
chiral superfield S as the dominant degree of freedom of low energy dynamics. Making a more 
general ansatz also containing gluon-gluon composites leads to a generalization and to two 
mixed supermultiplets in the low energy spectrum 0. Other generalizations are conceivable. 
This shows that the description of low energy dynamics in N=1 SYM theory by chiral effective 
actions is less rigorous than, for instance, in QCD where we know from the Goldstone theorem 
that the light states are the members of the pseudoscalar meson multiplet. 

A general consequence of spontaneous chiral symmetry breaking is the existence of a first 
order phase transition at nrig = 0. At this point the different ground states in ([17]) are degenerate 
and a coexistence of the corresponding phases is possible. In a mixed phase situation, as usual 
at first order phase transitions, the different phases are separated by “bubble wall” interfaces. 
The interface tension of the walls can be exactly derived from the central extension of the 
N = 1 SUSY algebra B The result is that the energy density of the interface wall is related 
to the jump of the gluino condensate by 


N 

e= 8^ l(AA)l “ (AA)a| • (25) 
Combining this with eq. (|T7D implies that the dimensionless ratio e/|(AA)| is predicted inde¬ 
pendently of the renormalization scheme. 

For transforming ([[7]) and ( |25| ) to lattice units we need, in fact, the value of aA lat at the 
particular values of interest of the lattice bare parameters (a is, as usual, the lattice spacing). 
Before performing the lattice simulations this is, of course, not known. An order of magnitude 
estimate can be obtained from pure gauge theory (at infinitely heavy gluino) by noting that 
both for N c = 2 and N c = 3 we have for the lowest glueball mass M [^] 

aA^r = ^ ■ (26) 

Assuming this approximate relation also at the critical line for zero gluino mass, we can use 
eqs. m-m for estimating orders of magnitudes. In the region where aM = 0(1) we obtain 


a 3 |(AA)i — (AA) 2 | = 0(1) , a 3 e = 0( 10" 1 ) . (27) 
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As these numbers show, the predicted first order phase transition is, in fact, strong enough for 
a relatively easy observation in lattice simulations. 

A basic question about supersymmetry is the way it is realized in quantum field theory. 
The phenomenology of broken supersymmetry is based on the assumption that in the super- 
symmetric extension of the Standard Model supersymmetry is not anomalous: there is only 
“soft breaking” due to mass terms. Since SUSY is broken by the known regularizations (for a 
discussion on this see |24|]), the question of a possible anomaly is legitimate. Non-perturbative 
effects may imply a supersymmetry anomaly Investigating the behaviour of SUSY Ward- 
Takahashi identities in the continuum limit of lattice regularization can help to give an answer 
to this question. 


2 Lattice formulation 

In order to define the path integral for a Yang Mills theory with Majorana fermions in the 
adjoint representation, let us first consider the familiar case of Dirac fermions [^(J. (Many 
of the notation conventions used in this paper are the same as in this reference.) Let us 
denote the Grassmanian fermion fields in the adjoint representation by -0“ and ip a x . Here Dirac 
spinor indices are omitted for simplicity and a stands for the adjoint representation index 
(a — 1,.., N% — 1 for SU(Y C ) ). The fermionic part of the Wilson lattice action is 

= te + ,wi + 7 j^+r^(i-7„«J } ■ (28) 

X fi=l 

Here K is the hopping parameter, the Wilson parameter removing the fermion doublers in 
the continuum limit is fixed to r = 1 and the matrix for the gauge-field link in the adjoint 
representation V x/J is defined from the fundamental link variables U XIX according to 

= Uv„[£7] = 2Tr(UlT a U v T b ) = . (29) 

The generators T a = |A a satisfy the usual normalization Tr (A a Ab) = 1. In the simplest case 
of SU(2) (N c = 2) we have, of course, T a = \r a with the isospin Pauli-matrices r a . The 
normalization of the fermion fields in (|28| ) is the usual one for numerical simulations. The full 
lattice action is the sum of the pure gauge part and fermionic part: 

S = S g + Sf. (30) 

The standard Wilson action for the SU(Y C ) gauge field S g is a sum over the plaquettes 

S g = P E( 1- IwReTr^) , (31) 

pi V 7V c ' 

with the bare gauge coupling given by (3 = 2 N c /g 2 . 






Ill order to obtain the lattice formulation of a theory with Majorana fermions let us note 
that out of a Dirac fermion field it is possible to construct two Majorana fields: 

A (1) = + Clf) , a< 2 > = ‘ (-V + cf ) 

with the charge conjugation matrix C. These satisfy the Majorana condition 

A W = A U)T C 


(32) 


0 = 1,2) . 


The inverse relation of (|32|) is 


ip = ^(A (1 > + iA< 2 >), y = c/ = ^(A(» - a< 2 >). 


(33) 


(34) 


In terms of the two Majorana fields the fermion action Sf in eq. (28) can be written as 

2 _ 4 


Sf = 5 E E { Af°AW“ - K y a““ y,, w (l + 7 „)A<’> + a“VJ, w (1 - 7„)a“‘ } . (35) 


X j =1 M=1 L 

For later purposes it is convenient to introduce the fermion matrix 

4 

, _ n , rrn = a a . _zy ’ 


Qyd,xc [u] = fiyx^dc y,x+p,( 1 + 7 /i)Vdc,xii + ^y+/i,x(l 7 r)^dc,yy, 


M= 1 


(36) 


Here, as usual, /z denotes the unit vector in direction /i. In terms of Q we have 

1 


7= E 5 E E v w • 


1 — (j)rf 

xc,yd ^ j =1 xc,yd 

and the fermionic path integral can be written as 

2 


(37) 


e~ Sf = 


e -yw = det Q = n 

i=i 




(38) 


This shows that the path integral over the Dirac fermion is the square of the path integral over 
the Majorana fermion and therefore 


- 5 aqa = ± D det g 


(39) 


As one can see here, for Majorana fields the path integral involves only [eiA^] because of the 
Majorana condition in (|33|). 

The relation ([59]) leaves the sign on the right hand side undetermined. A unique definition 
of the path integral over a Majoran fermion held, including the sign, is given by 


e -iAQA = 


e -iAMA 


= Pf (M) 


where M is the antisymmetric matrix defined as 

M = CQ — -M t 


(40) 


(41) 
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The square root of the determinant in eq. fl39|) is a Pfaffian. This can be defined for 
a general complex antisymmetric matrix M a p = —Mp a with an even number of dimensions 
(1 < a, (3 < 2N ) by a Grassmann integral as 


Pf (M) = 


e _ —— —e ai01 ... aN !3 N M a ^ ... M aNl 3 N 


(42) 


Here, of course, [d(j)\ = d(p 2 N ■ ■ ■ d<f> i, and e is the totally antisymmetric unit tensor. 

It is now clear that the fermion action for a Majorana fermion in the adjoint representation 
A“ can be defined by 

4 - - 

}. (43) 


Sf = xA«A = W{ A“A“ - K y [X a , +ll V al „(l + 7 „)A‘ + X x V^(l - 7 PA ‘ +(i 

Z ^ X fj ,=1 


This together with (^Dj)-(^Tj) gives a lattice action for the gauge theory of Majorana fermion 
in the adjoint representation. In order to achieve supersymmetry one has to tune the hopping 
parameter (bare mass parameter) K to the critical value K cr (/3) in such a way that the mass 
of the fermion becomes zero. 

The path integral over A is defined by the Pfaffian Pf(CQ) = Pf(Af). By this definition the 
sign on the right hand side of eq. (]39|) is uniquely determined. The determinant det(Q) is real 
because the fermion matrix in (|36|) satisfies 


Q ] = IsQlz , 


Q = 75 Q = Q ] • 


(44) 


Moreover one can prove that det(Q) = det(Q) is always non-negative. This follows from the 
relations 

CQC - 1 = Q t , BQB- 1 = Q t , (45) 

with the charge conjugation matrix C and B = C 75 . It follows that every eigenvalue of Q and 
Q is (at least) doubly degenerate. Therefore, with the real eigenvalues A* of the Hermitean 
fermion matrix Q , we have 

det(Q) = det(Q) = - 0- (46) 


Since according to the above discussion 


det(Q) = det(M) = [Pf(A /)] 2 , 


(47) 


the Pfaffian Pf(ilf) has to be real - but it can have any sign. 

The relation between Majorana- and Dirac-fermions can also be exploited for the calculation 
of expectation values of Majorana fermion fields. For the Dirac fermion fields i/j, ifj we have, as 
is well known, 


• • -VvAJ = Z 


7-1 




= z 


-1 


e ~s g [u] det Q[U] 


c ZiZ2---Z n 

~ yi y2---yn' 




1-1 


(48) 
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with the antisymmetrizing unit tensor e and 


Z = 


i-Sg[U]-ipQip _ 


e -S,P] det Q[t7] . 


(49) 


In order to express the expectation value of Majorana fermion fields by the matrix elements 
of the propagator one can use again the doubling trick: one can identify, for instance, 

A = A (1) and introduce another Majorana field A*- 2 ), in order to obtain a Dirac held. Then from 
eqs. (|32|) and ([48]) we obtain 


^y2 ^x 2 ' ' ' ^Xvt / Z 


IJn /x Xn 


7 ~m I[dU] e~ Sa[u] Pf(M[D])(detg[f/])- 1 




(V’w + C ^ yi )(^ xi (VVn + ^ x n C ) 2 


where now 


Z M = j{dU}e~ Sa[u] Pf [M[U]) . 
The simplest example is n — 1 when we have with Q = Q[U} 


(50) 

(51) 


^y^x) Z m 


= Z 


-1 

M 


e~ Salu] Pf (M[U]) \ {Qyt + C- l Q~ l y C} 


e- s A u } pf (M[U\) Q~l . 


(52) 


The important case n = 2 can be expressed by six terms: 


K\ X2 \ yi \ Xl ) = Zl} J[ciU]e~ s Pf {M[U}) 


i22/n-i /n-i i f)- 1 n _l r* -1 ^ 2122 (l~ l C 1 

\ C yiV2^V Z 2 X 2 ^ C ZlJ/2^* Ziyi^ Z2X 2 ^yi ' KJ X 2 c y lX2 ^ZlXl Ka 6z 2 y2^'! 

2122 


U2 


n~ l n~ 1 c c - o~ 1 c — n~ l c 1 (cmt 

T ^n U I2 C IlI2 l 1i2 1 !) 1 'i:z2!)2 W !/2 X! C yiXl^Z 1X 2^Z2y2 V2 ^ X2 y2X2 ^ 21*1 ^z 2 yi '-'yif ■ \ 00 J 

The indices on the charge conjugation matrix C show how the Dirac indices have to be con¬ 
tracted. 

Another way of expressing the expectation values of Majorana fermions is to consider the 
path integral with an external source J x 


Z M [J] = J [dU}e~ Salu] 

After a shift of the Grassmanian integration variable 

A' = A + JM~ l 

one obtains 

Z M [J\ = I{dU] e - s ‘M j[d\']e-i yMX '-i JM - 1J 


e —±AMA—JA 


e~ Sa[u] Pf (M) e" 2 


-ijM-D 


(54) 


(55) 


(56) 


Expectation values of the form (|52|) , (|53|) can now be obtained by differentiation with respect 
to the source J. 
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2.1 Improved actions 


The lattice action defined by ( pOj ), ( pll) and (J43P is not unique. One expects that there is a wide 
class of actions belonging to the same universality class which all reproduce the same theory 
in the continuum limit. This ambiguity can be used for improvement of some desired property, 
for instance, faster approach to the continuum limit and/or less symmetry breaking. 

Direct improvement of the supersymmetry of the lattice action seems to be rather hard 
in gauge theories because the gauge field and the fermions are introduced on the lattice very 
differently. Nevertheless, in Wess-Zumino type scalar-fermion models it is possible to achieve 


perfect supersymmetry on the lattice with respect to supersymmetry type transformations |28 
(See also 


Another possibility of improvement is to apply a lattice formulation for fermions where the 
zero gaugino mass is achieved without fine tuning. The preferred way to do this is based on 
domain wall fermions |30| where an auxiliary fifth dimension is introduced which contains a 
four dimensional domain wall supporting a massless fermion. A similar alternative is based on 
overlap fermions [[H], |32| . In these formulations an additional advantage is the improved chiral 
symmetry. In fact, these lattice fermions are prominent representatives of the Ginsparg-Wilson 
realization of chiral symmetry on the lattice . 

An important advantage of Ginsparg-Wilson fermions is due to their clean spectral proper¬ 
ties: for zero mass the eigenvalues are in the complex plane on a circle crossing the real axis 
at zero. For non-zero mass the crossing is shifted by the mass in lattice units. The product of 
eigenvalues (determinant or Pfaffian) cannot be negative if the mass is non-negative. This is in 
contrast with Wilson-type fermions where the fluctuation of eigenvalues implies the possibility 
of negative eigenvalues also for (small) positive masses. This is the origin of a possible “sign 
problem” with Wilson-type lattice fermions (see section |3.2|) . 

Most of the numerical simulations up to now used the simple Wilson-type lattice action in 
(|43|) . These will be discussed in detail in sections fif|| Recently there has been, however, also a 
Monte Carlo study of SYM theory based on the domain wall fermion action |34j . In this work 
an evidence was found for the formation of a gaugino condensate in the chiral limit at zero 
gaugino mass. 


3 Algorithm for numerical simulations 

In order to perform Monte Carlo simulations of SYM theory one needs a positive measure on 
the gauge field which allows for importance sampling of the path integral. Therefore the sign 
of the Pfaffian can only be taken into account by reweighting. According to ([|7]) the absolute 
value of the Pfaffian is the non-negative square root of the determinant therefore the effective 
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gauge field action is [35]: 


S CV = PJ2^ 

pi 


-Tr a 


pi 


log det Q[U) 


(57) 


The factor \ in front of log det Q shows that we effectively have a flavour number Nf = A of 
adjoint fermions. The omitted sign of the Pfaffian can be taken into account by reweighting 
the expectation values according to 


( 0 ) = 


(O signPf (M))cv 
(signPf(M)) cy 


(58) 


where (.. .)cv denotes expectation values with respect to the effective gauge action Scv • This 
may give rise to a sign problem which will be discussed in section |3.2| . 

The fractional power of the determinant corresponding to (|)7|) can be reproduced, for in¬ 
stance, by the hybrid molecular dynamics algorithm [[36[] which is, however, a finite step size 
algorithm where the step size has to be extrapolated to zero. This adds another necessary 
extrapolation to three: to infinite volume, to small gaugino mass and to the continuum limit. 
An “exact” algorithm where the step size extrapolation is absent is the two-step multi-boson 
(TSMB) algorithm |37| . 


3.1 Two-step multi-boson algorithm 


The multi-boson algorithms for dynamical (“unquenched”) fermion simulations |3£| are based 
on polynomial approximations as 


|det(Q)| iV/ = |det(Q t Q)| 


Nf/ 2 


(59) 


detP n (QtQ) ’ 

where Nf is the number of flavours (in our case Nf = |) and the polynomial P n satisfies 

lim P n (x) = x ~ Nf / 2 (60) 

n— kx) 

in an interval [e, A] covering the spectrum of Q'Q = Q 2 . (Here Q is the hermitean fermion 
matrix defined in (|44|).) 

For the multi-boson representation of the determinant one uses the roots of the polynomial 
rj, U = l,...,n) 

n 

PniQ'Q) = Pn(Q 2 ) = r 0 n (Q 2 - rj) . (61) 

Assuming that the roots occur in complex conjugate pairs, one can introduce the equivalent 
forms 


Pn(Q 2 ) = 'M ]['(Q ± h ) 2 + V 2 ] = r 0 P(Q - p’)(Q - Pi) 


(62) 


3 =1 


3 = 1 
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where r 3 = (jij + ivj) 2 and p 3 = /ij + iv 3 . With the help of complex boson (pseudofermion) 
fields one can write 


n r \ n 

n det [(Q - Pj)(Q - Pj)}~ 1 tx / [d$] exp ^ ^ ^ $+ [(Q - p*){Q - Pj)\ yx § jx \ . (63) 

j =i l 3=t x y I 


Since for a finite polynomial of order n the approximation in (| 6 Dj) is not exact, one has to 
extrapolate the results to n —> oo. 

The difficulty for small fermion masses in large physical volumes is that the condition number 
A/e becomes very large (10 4 — 10 6 ) and very high orders n = 0(1O 3 ) are needed for a good 
approximation. This requires large storage and the autocorrelation of the gauge configurations 
becomes very bad since the autocorrelation length is a fast increasing function of n. One can 
achieve substantial improvements on both these problems by introducing a two-step polynomial 
approximation: 

P $(. x ) P n 2 { x ) = x ~ Nf/2 , x e [e, A] . (64) 

The multi-boson representation is only used for the first polynomial P// which provides a first 
crude approximation and hence the order ri\ can remain relatively low. The correction factor 
P/T j s realized in a stochastic noisy correction step with a global accept-reject condition during 
the updating process. In order to obtain an exact algorithm one has to consider in this case 
the limit n 2 —■► oo. 

In the two-step approximation scheme for Nf flavours of fermions the absolute value of the 
determinant is represented as 


|det(Q)| V/ ~ 


1 

det Pni (Q 2 ) det PnJ (Q 2 ) 


(65) 


The multi-boson updating with n\ scalar pseudofermion fields is performed by heatbath and 
overrelaxation sweeps for the boson fields and Metropolis (or heatbath and overrelaxation) 
sweeps over the gauge field. After an update sweep over the gauge field a global accept-reject 
step is introduced in order to reach the distribution of gauge field variables [U] corresponding 
to the right hand side of (^5|). The idea of the noisy correction is to generate a random vector 
77 according to the normalized Gaussian distribution 


e -r,'P™(Q[U] 2 )r, 

J[drj\e-^ p ^^ u l 2 ) 77 ’ 


( 66 ) 


and to accept the change [U] —> [U'\ with probability 


min {1,34(77; P'\ <- [U])} , 


(67) 


where 

\U'\ <- [£/]) = exp{-,tp«(Q[C/f), + P£(Q[Vf)r,} . ( 68 ) 
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The Gaussian noise vector r/ can be obtained from ;/ distributed according to the simple 
Gaussian distribution 




—77't 77 ' 


by setting it equal to 


SW]e 


1 = p^mufrh 1 . 


(69) 


(70) 


In order to obtain the inverse square root on the right hand side of ([71]), one can proceed 
with a polynomial approximation 


P$( x ) - Pm( x ) b x 6 [0,A] . 


(71) 


This is a relatively easy approximation because P^{x)~% is not singular at x — 0, in contrast 
to the function x~ N ^ 2 . A practical way to obtain P ® is to use some approximation scheme 
for the inverse square root. The best possibility is to use a Newton iteration 


Pm = i ( Pi?' + 


( 3 ) 


p( 3) p(2) 


k = 0,1,2,... 


(72) 


The second term on the right hand side can be evaluated by a polynomial approximation as 
for P in ( [64]) with Nf = 0 and P W —»■ Pjf' 1 P( 2 \ The iteration in (|72[) is fast converging and 

/o\ 

allows for an iterative procedure stopped by a prescribed precision. A starting polynomial Pq 
can be obtained, for instance, from (| 64 |) with Nf —> — \Nf or by the formula |39| 


1 K 
p(3) ~ _ V 

K ^ 


8=1 


P( 2) 


COS 


2 7T 

2 K 


:s-l)+sm 2 ^(s-l) 


(73) 


The TSMB algorithm becomes exact only in the limit of infinitely high polynomial order: 
n 2 —■► 00 in (^j4[). Instead of investigating the dependence of expectation values on n 2 by per¬ 
forming several simulations, it is better to fix some high order n 2 for the simulation and perform 
another correction in the “measurement” of expectation values by still finer polynomials. This 
is done by reweighting the configurations. 

Before discussing this measurement correction for general number of flavours Nf let us note 
that Nf = 2, and any other even number, is a special case. This is because, due to P^^x) ~ 
[ x p( 1 )( a; )]- 1 , the noisy correction can be performed in this case by an iterative inversion, as 
first introduced in [401. Consequently, one can proceed without a measurement correction. In 
another multi-boson scheme JIT]] (polynomial hybrid Monte Carlo) a measurement correction 
is needed but for even Nf it can be performed by an iterative inversion. 

The measurement correction for general Nf is based on a polynomial approximation P^ 
which satisfies 

lim \ x )PnX x )PnXx) = x ~ Nf/2 , X G [e', A] . (74) 


714—XX) 


The interval [e 7 , A] can be chosen by convenience, for instance, such that e' = 0, A = A max , 
where X max is an absolute upper bound of the eigenvalues of Q^Q = Q 2 . (In practice, instead 
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of e' — 0 , it is more effective to take e' > 0 and determine the eigenvalues below e' and the 
corresponding correction factors explicitly |4^].) In this case the limit 714 —> 00 is exact on an 
arbitrary gauge configuration. For the evaluation of one can use n 4 -independent recursive 
relations (43j, which can be stopped by observing the required precision of the result. After 
reweighting the expectation value of a quantity O is given by 

(Oexp{ v tll-P! l ?(Q<Q ' l } v }) u , 1 


(O) = 


(exp {i;t[l - PS$(QHQ)]q})u r 


(75) 


where // is a simple Gaussian noise like 7/ in (p9|) . Here (.. .)u,ri denotes an expectation value 
on the gauge field sequence, which is obtained in the two-step process described in the previous 
subsection, and on a sequence of independent 7 /’s. The expectation value with respect to the ■//- 
sequence can be considered as a Monte Carlo updating process with the trivial action S v = 77 b/. 
The length of the //-sequence on a fixed gauge configuration can, in principle, be arbitrarily 
chosen. In praxis it has to be optimized for obtaining the smallest possible errors. If the 
second polynomial gives a good approximation the correction factors do not practically change 
the expectation values. In this reweighting step the sign of the Pfaffian can also be included 
according to (|58|) and then one has 


(0) = 


(O signPf(M) exp {?/t[l - P${Q*Q)]ri})u,- 
(signPf(M) exp {r/t[l - P,u (Q^Q)]v}) 


(76) 


A useful generalization of the two-step multi-boson algorithm is to allow for the fermion 
action in the multi-boson update step to differ from the true fermion action one wants to 
simulate [|44|. For instance, for improved actions one can take a simplified version in the multi¬ 
boson update step. Random variation of the parameters in the update step may also be used 
to decrease autocorrelations. 


3.2 The “sign problem” 

The Pfaffian resulting from the Grassmannian path integrals for Majorana fermions (]4^) is an 
object similar to a determinant but less often used. As shown by ([42]), Pf(M) is a polynomial of 
the matrix elements of the 2iV-dimensional antisymmetric matrix M = —M T . Basic relations 
are |45[ 

M = P T JP , Pf(Af) = det(P) , (77) 

where J is a block-diagonal matrix containing on the diagonal 2 ® 2 blocks equal to e = ia 2 


and otherwise zeros. Let us note that from these relations the second equality in eq. 
immediately follows. 

The form of M required in ( |77| ) can be achieved by a procedure analogous to the Gram- 
Schmidt orthogonalization and, by construction, P is a triangular matrix (see ii). This gives 
a numerical procedure for the computation of P and the determinant of P gives, according to 
(|77|), the Pfaffian Pf(M). Since P is triangular, the calculation of det(P) is, of course, trivial. 
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This procedure can be used for a numerical determination of the Pfaffian 
|46| . For an example see figure [L] which also shows a sign change of Pf(M) as 
hopping parameter K. 


on small lattices 
a function of the 


lattice: 4040408 



■^valence 


Figure 1: The absolute value of the pfaffian on a 4 3 -8 configuration 
as a function of the hopping parameter in M (and Q). Open points 
stand for Pf(M) > 0, full ones for Pf(M) < 0. 


On lattices larger than, say, 4 3 • 8 the computation becomes cumbersome due to the large 
storage requirements. This is because one has to store a full 0 0 0 matrix, with O being the 
number of lattice points multiplied by the number of spinor-colour indices (equal to 4(Ay 2 —1) for 
the adjoint representation of SU(1V C )). The difficulty of computation is similar to a computation 
of the determinant of Q with L£/-decomposition. 

Fortunately, in order to obtain the sign of the Pfaffian occurring in the measurement 
reweighting formulas (|58|), ([76]) one can proceed without a full calculation of the value of 
the Pfaffian. The method is to monitor the sign changes of Pf(M) as a function of the hopping 
parameter K. According to (|46|) , the hermitean fermion matrix for the gaugino Q has doubly 
degenerate real eigenvalues therefore 

n/2 

det M = det Q = , (78) 

i= 1 


where A i denotes the eigenvalues of Q. This implies 


0/2 

iPf(M)i = n iA f 

i =1 


Q/2 

pf(M) = n Ai. 

i=1 


(79) 


The first equality trivially follows from (]47|). The second one is the consequence of the fact 
that Pf(M) is a polynomial in K which cannot have discontinuities in any of its derivatives. 
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Therefore if, as a function of K , an eigenvalue A* (or any odd number of them) changes sign 
the sign of Pf(M) has to change, too. Since at K = 0 we have Pf(M) = 1, the number of sign 
changes between K = 0 and the actual value of K, where the dynamical fermion simulation is 
performed, determines the sign of Pf(M). This means that one has to determined the flow of 
the eigenvalues of Q through zero [[IT]. Examples of the spectral flow taken from data of the 
Monte Carlo simulations of the DESY-Munster Collaboration jl2] are shown in figure ^]. 
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Figure 2: The spectral flow of the hermitean fermion matrix Q for some specific 
configurations on 6 3 • 12 lattice at (3 = 2.3. The value of K in the simulation is 
displayed by a vertical line. 


The spectral flow method is well suited for the calculation of the sign of the Pfaffian in SYM 
theory. An important question is the frequency and the effects of configurations with nega¬ 
tive sign. A strongly fluctuating Pfaffian sign is a potential danger for the effectiveness of the 
Monte Carlo simulation because cancellations can occur resulting in an unacceptable increase 
of statistical errors. The experience of the DESY-Munster Collaboration shows, however, that 
below the critical line K cr (f3 ) corresponding to zero gaugino mass (mg = 0) negative Pfaffians 
practically never appear ||2], [48], [49] . Above the critical line several configurations with negative 
Pfaffian have been observed but their role has not yet been cleared up to now. Since supersym¬ 
metry is expected to be realized in the continuum limit at mg = 0, the negative signs of the 
Pfaffian can be avoided if one takes the zero gaugino mass limit from rrig > 0 corresponding to 
K < K cr . In this sense there is no “sign problem” in SYM with Wilson fermions which would 
prevent a Monte Carlo investigation. 
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The presence or absence of negative Pfaffians in a sample of gauge configurations produced 
in Monte Carlo simulations can be easily seen even without the application of the spectral flow 
method. In case of sign changes the distribution of the reweighting factors in (f75[) shows a 
pronounced tail reaching down to zero |50|. If this tail is absent, as in the example of a run 


fTj| shown by figure |Sj. then there are no negative Pfaffians. Alternatively, one has to observe 



0 0.25 0.5 0.75 1 1.25 1.5 1.75 

Figure 3: The distribution of reweighting factors at (3 = 2.3, K = 
0.1940 on 12 3 ■ 24 lattice. 


the distribution of smallest eigenvalues of Q 2 which also shows a tail down to zero if there are 
negative Pfaffians in the Monte Carlo sample. The absence of a tail shows that there are no 
negative Pfaffians. This is illustrated by figure ^ belonging to the same configuration sample 
as figure ||. 

Concerning this “sign problem” let us note that a very similar phenomenon exists also in 
QCD because the Wilson-Dirac determinant of a single quark flavour can also have a negative 
sign. Under certain circumstances the sign of the quark determinant plays an important role. 
This is the case, for instance, at large quark chemical potential in a QCD-likc model with SU(2) 
colour and staggered quarks in the adjoint representation which has recently been studied by the 
DESY-Swansea Collaboration []5U|. This investigation also revealed an interesting feature of the 
TSMB algorithm, namely its ability to easily change the sign of eigenvalues of the hermitean 
fermion matrix (and hence the sign of the determinant or Pfaffian). This is in contrast to 
algorithms based on finite difference molecular dynamics equations as, for instance, the hybrid 
Monte Carlo |5lJ or HMD |36| algorithms. 
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Figure 4: The distribution of smallest eigenvalues of Q 2 at (3 = 
2.3, K = 0.1940 on 12 3 • 24 lattice. 


4 Discrete chiral symmetry breaking 


It has been discussed in section [TTl] that in SYM theory with gauge group SU(N C ) the discrete 
global chiral symmetry Z 2 n c is expected to be spontaneously broken by a non-zero gaugino 
condensate (XX) ^ 0. The consequence of this chiral symmetry breaking pattern is the existence 
of a first order phase transition at zero gaugino mass ntg = 0. 

In case of the lattice formulation based on the lattice action in (|30|) - (l3l|) , ([h|) the bare 
parameters are the gauge coupling (3 = 2 N c /g 2 and the hopping parameter K. In the plane 
(/3, K) there is a line corresponding to zero gaugino mass K cr (/3 ) where a first order phase 
transition occurs. Therefore the expected phase structure is the one shown in figure ^|. In 
fact, the theoretical expectation of a genuine phase transition refers to the continuum limit 
(3 —> oo. It is possible that for finite (3 there is a cross-over which becomes a real first order 
phase transition only at (3 = oo. 

The DESY-Miinster Collaboration performed a first lattice investigation of gaugino conden¬ 
sation in SYM theory with SU(2) gauge group ]5!J. The distribution of (AA) has been studied 
for fixed gauge coupling (3 = 2.3 as a function of the hopping parameter K, which determines 
the bare gaugino mass, on 6 3 ■ 12 lattice. 

A Erst order phase transition shows up on small to moderately large lattices as metastability 
expressed by a two-peak structure in the distribution of the gaugino condensate. By tuning 
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Figure 5: Expected phase structure of Yang-Mills theory with a Ma- 
jorana fermion in adjoint representation in the (f3, K)-plane. The 
dashed-dotted line K = K cr (/3) is a first order phase transition (or 
cross-over) at zero gaugino mass. 


K one should be able to achieve that the two peaks are equal (in height or area). This is a 
possible definition of the phase transition point in finite volumes. By increasing the volume 
the tunneling between the two ground states becomes less and less probable and at some point 
practically impossible. 

The observed distributions are shown in figure |6j. One can see that the distributions cannot 
always be described by a single Gaussian which would correspond to a single phase. Parameters 
of possible fits are collected in table |l|. As the table shows, in the region 0.195 < K < 0.1975 
only two-Gaussian fits work. Outside this region single Gaussian fits describe the data well. 
For increasing K (decreasing bare gaugino mass) the weights shift from the Gaussian at larger 
p = (XX) to the one with smaller p, as expected. The two Gaussians represent the contributions 
of the two phases on this lattice. The position of the phase transition (or cross-over) on the 
6 3 • 12 lattice is at K 0 = 0.1955 ± 0.0005. According to table |l] the jump of the order parameter 
is A p = p 1 — p 2 — 0.15. 


The two-phase structure can also be searched for in pure gauge field variables as the pla- 
quette or longer Wilson loops. It turns out that the distributions of Wilson loops are rather 
insensitive. They can be well described by single Gaussians with almost constant variance in 


the whole range 0.19 < K < 0.2 [52 
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Figure 6: The probability distributions of the gaugino condensate p = (XX) for 
different hopping parameters at /3 = 2.3 on 6 3 • 12 lattice. The dased lines show 
the Gaussian components. 
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Tabic 1: Fit parameters of the distributions of gaugino condensate corresponding 
to figure |||. The weight of the first peak is pi, the mean values are /j^ 2 and the 
common width of the Gaussians is a. The statistical errors in last digits are given 
in parentheses. 


K 

Pi 

Pi 

p2 

(7 

X 2 / d.o.f. 

0.19 

1.0 

11.0023(26) 

- 

0.0423(16) 

27.9/20 

0.1925 

1.0 

10.8807(30) 

- 

0.0524(17) 

25.9/20 

0.195 

0.89(7) 

10.762(30) 

10.608(30) 

0.066(7) 

16.5/18 

0.196 

0.35(7) 

10.722(11) 

10.588(11) 

0.073(3) 

5.7/18 

0.1975 

0.26(5) 

10.626(17) 

10.484(17) 

0.056(4) 

19.5/18 

0.2 

0.0 

- 

10.3363(37) 

0.0562(18) 

21.4/20 


As shown by figure and table [I|, on this small lattice the two peaks corresponding to two 
phases are not well separated. It is possible that at this (3 value there is no genuine first order 
phase transition at all - only a cross-over. The question of first order phase transition versus 
cross-over can be decided by studying the volume dependence on larger lattices. After extracting 
the jump of the gaugino condensate in the large (infinite) volume limit one has to study its 
behaviour in the continuum limit (3 —> oo along the line K cr {f3 ) in figure || The renormalized 
gaugino condensate can be obtained in the continuum limit by additive and multiplicative 
renormalizations: 


(X x X x )R(p) Z(ap) (X x X x ) b 0 (an) 


(80) 


The presence of the additive shift in the gaugino condensate b 0 (ap) implies that the value of 
its jump at ntg = 0 is easier available than the value itself. 

The renormalization factor Z in (^U|) is expected to be of order 0(1) at bare parameter values 
which can be reached by numerical simulations. For obtaining it one can use non-perturbative 
renormalization schemes (for a recent review see |[!|). Performing perturbative calculations in 
lattice regularization can also give useful information 


5 Confinement and particle spectrum 

5.1 Fundamental static potential 

The potential between static colour sources in gauge field theory is a physically interesting 
quantity because it is characteristic for the dynamics of the gauge field. The external static 
soursec can be put in any representation of the gauge group. If the sources are in the fun¬ 
damental representation we have to do with the fundamental static potential between static 
quarks. 
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For a model containing dynamical matter fields in the fundamental representation, as is the 
case for QCD, the charge of static quarks will be screened. The potential then approaches a 
constant at large distances. The string tension a, which is the asymptotic slope of the potential 
for large distances, vanishes accordingly. On the other hand, if only matter fields in the adjoint 
representation of the gauge group are present, as in the case of N=1 SYM theory, confinement 
of static quarks and a positive a are expected. 

The DESY-Miinster Collaboration has determined the static quark potential and the string 
tension for N=1 SUSY Yang-Mills theory with gauge group SU(2) in Monte Carlo simulations 
[|42f . The starting point of numerical work are expectation values of rectangular Wilson loops 
(W(R, T)) of spatial length R and time length T. From the Wilson loops the potential can be 
found via 

V(R) = lim V(R,T ) , (81) 

1 —>oo 

where 

V(R, T) = log (W(R, T)) - log (W(R, T + 1)) . (82) 

The potential V(R) is obtained through a fit of the form 

V(R, T) = V(R) + ci (R)e~ C2(R)T . (83) 

The string tension a is finally obtained by fitting the potential according to 

Oi 

V(R) = V 0 -- + aR. (84) 

it 

An example of the static quark potential on 12 3 • 24 lattice at j3 — 2.3, K = 0.1925 is shown 
in figure |7]. A collection of the results for the string tension on lattices with spatial extension 
L is: 


ay/a = 

0.22(1) 

for K = 0.1900, L = 8, 


ay/a = 

0.21(1) 

for K = 0.1925, L = 8, 


ay/a = 

0.17(1) 

for K = 0.1925, L = 12. 

(85) 


The string tension in lattice units is decreasing when the critical line is approached, as it should 
be. This is mainly caused by the renormalization of the gauge coupling due to virtual gluino 
loop effects which are manifested by decreasing lattice spacing a. From a comparison of the 
L — 8 and L — 12 results one sees that finite size effects still appear to be sizable. This has to 
be expected because we have for the spatial lattice extension L = 12 a the result Ly/a ~ 2.1. 
In QCD with /a ~ 0.45 GeV this would correspond to L ~ 1 fm. Although we are dealing 
with a different theory where finite size effects as a function of Ly/a are different, for a first 
orientation this estimate should be good enough. 

For the ratio of the scalar glueball mass m(0 + ), to be discussed below, and the square root 
of the string tension the results are: 

m(0 + )/V° = 3.4(7) for K = 0.1900, L = 8, 
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Lattice: 12 3 x24, K = 0.1925 



Figure 7: The static quark potential V(R) on a 12 3 • 24 lattice at 
(3 = 2.3, K = 0.1925. The line is a fit with a Coulomb plus a linear 
term, fitted over the range 1 < R < 6. 

m(0 + )/y/a = 3.0(4) for K = 0.1925, L = 8, 

m{ 0 + )/v^ = 3.1(7) for K = 0.1925, L = 12. (86) 

The uncertainties are not very small, but the numbers are consistent with a constant indepen¬ 
dent of K in this range. 

5.2 Supersymmetry multiplets? 

The non-vanishing string tension implies that the Yang-Mills theory with gluinos is confining. 
Therefore the asymptotic states are colour singlets, similarly to hadrons in QCD. The structure 
of the light hadron spectrum is closest to the (theoretical) case of QCD with a single flavour of 
quarks where the chiral symmetry is broken by the anomaly. 

Since both gluons and gluinos transform according to the adjoint (for SU(2) gauge group 
triplet) representation of the colour group, one can construct colour singlet interpolating fields 
from any number of gluons and gluinos if their total number is at least two. Experience in QCD 
suggests that the lightest states can be well represented by interpolating fields built out of a 
small number of constituents. Simple examples are the glueballs known from pure Yang-Mills 
theory and gluinoballs which are composite states made out of gluinos. Mixed gluino-glueball 
states can be composed of any number of gluons and any number of gluinos, in the simplest 
case just one of both. 

In general, one has to keep in mind that the classification of states by some interpolating 
fields has only a limited validity, because this is a strongly interacting theory where many 
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interpolating fields can have important projections on the same state. Taking just the simplest 
colour singlets can, however, give a good qualitative description. 

In the supersymmetric limit at zero gluino mass rrig = 0 the hadronic states should occur 
in supermultiplets. This restricts the choice of simple interpolating field combinations and 
leads to low energy effective actions in terms of them [|T^, |2]J. For non-zero gluino mass the 
supersymmetry is softly broken and the hadron masses are supposed to be analytic functions 
of rrig. The linear terms of a Taylor expansion in rrig are often determined by the symmetries 
of the low energy effective actions m- 

The effective action of Veneziano and Yankielowicz [18] in eq. (23) describes a chiral super- 
multiplet consisting of the 0“ gluinoball fj, the 0 + gluinoball jo and a spin | gluino-glueball 
X■ There is, however, no a priori reason to assume that glueball states are heavier than the 
members of this supermultiplet. Therefore Farrar, Gabadadze and Schwetz [^T[ proposed an ef¬ 
fective action which includes an additional chiral supermultiplet. This consists of a 0 + glueball, 
a 0“ glueball and another gluino-glueball. The effective action allows mass mixing between the 
members of the two supermultiplets. 

The spectrum of SYM theory with SU(2) gauge group has been studied in numerical sim¬ 
ulations by the DESY-Mimster Collaboration |I2|, |5S], |57]]. The glueball states as well as the 
methods to compute their masses in numerical Monte Carlo simulations are well known from 
pure gauge theory (see [|23|). The obtained masses for the 0 + glueball in lattice units are 


am(0 + ) 

= 0.95(10) 

am( 0 + ) 

= 0.85(6) 

am(0 + ) 

= 0.75(6) 

am(0 + ) 

= 0.63(5) 

am(0 + ) 

= 0.53(10) 


for K = 0.1800, L = 6, 
for K = 0.1850, L = 6, 
for K = 0.1900, L = 8, 
for K = 0.1925, L = 8, 
for K = 0.1925, L = 12. 


(87) 


As before, L denotes the spatial extension of the lattice. (The time extension has always been 
T = 2 L.) 

In addition to the J p = 0 + glueball one can also search for the pseudoscalar 0 - glueball. 
The masses in lattice units turned out to be 


am(0~) 

= 1-5(3) 

for K = 0.1850, L = 6, 


am( 0“) 

= 1.45(10) 

for K = 0.1900, L = 6, 


am( 0 _ ) 

= 1-3(1) 

for K = 0.1925, L = 6, 


am( 0“) 

= 1-1(1) 

for K = 0.1925, L = 8. 

(88) 


The pseudoscalar glueball appears to be roughly twice as heavy as the scalar one. This is 
similar to pure SU(2) gauge theory, where m(0 _ )/m(0 + ) = 1.8(2) 

As discussed above, besides the glueballs made out of gluons one also has to consider 
gluinoballs. Examples of simple colourless composite fields can be constructed from two gluino 
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fields: A^A* and A^A^. These correspond to the gluinoball states mentioned above: fj and / 0 , 
respectively, which are described by the effective action (p3|). 

The masses of fj and / 0 can be extracted from correlation functions as 


r sj (x,y) = (TVsc {rg- 1 } - 2TV„{rQ- 1 r<?- 1 }), 


(89) 


where T G { 1 , 75 } and Q 1 is the glnino propagator. Note that the factor of two is originating 
from the Majorana character of the glninos. In analogy with the flavour singlet meson in QCD 


the correlator in fl89[) consists of a connected and a disconnected part. The disconnected part 
can be calculated using the volume source technique |58| . 


I 11 case of the /o particle the disconnected and the connected parts are of the same order of 
magnitude whereas fj is dominated by the connected part. The disconnected part has a much 
worse signal to noise ratio than the connected one. This leads to a larger error on the /o mass 
as compared to the f] mass. 

In the low energy effective action of Farrar, Gabadadze and Schwetz m there is a possible 
non-zero mixing between the states in two light supermultiplets. I 11 particular there can be 
mixing of the / 0 gluinoball and the 0 + glueball which have identical quantum numbers. The 
numerical simulations show that this mixing is small |4^, [57| . 

The low mass supermultiplets containing J p = 0 _ , 0 + states have to be completed by a spin 
A state. The corresponding composite fields can constructed from the glnino field and the field 
strength tensor of gluons. The simplest example is WA |57|: 


$ = Tr [F^ A] . 


(90) 


On the lattice one has to use, of course, an appropriate construction for F^ u out of the link 
variables U XfX . Composite fields with the same quantum numbers can also be built up from 
three glnino fields |56| , [57]. A summary of the results of the DESY-Munster Collaboration 
about the spectrum of light composite states is shown in figure ||. 

As the figure shows, the low mass states appear to build two groups which may correspond 

For a firm conclusion more numerical work is needed in larger 
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volumes and for several values of the gauge coupling f3 allowing for a continuum extrapolation. 


6 Supersymmetric Ward-Takahashi identities 

An important feature of lattice regularization is that some symmetries are broken for non-zero 
lattice spacing and are expected to be recovered in the continuum limit. The details of the 
lattice formulation, which may also influence the degree of symmetry breaking, are not relevant 
in the continuum limit because of the universality of critical points. A basic set of symmetries 
broken by the lattice and restored in the continuum limit is the (Euclidean) Lorentz symmetry 
including rotations and translations. It is clear that on any regular lattice these symmetries 
are always broken. Internal symmetries as, for instance, global chiral symmetry are sometimes 
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Figure 8: The lightest bound state masses in lattice units as func¬ 
tion of the bare gluino mass parameter 1/K for fixed (3 = 2.3. The 
shaded area at K = 0.1955(5) is where zero gluino mass and super- 
symmetry are expected. In this figure the notations of |5^, [5j/ 
are used: a-r/ = ?) and a-f 0 = / 0 . 


broken sometimes conserved on the lattice, depending on the actual formulation. From the 
point of view of symmetry realization supersymmetry is expected to behave similarly to the 
Lorentz symmetry: at finite lattice spacing it is broken but it becomes restored in the continuum 
limit. This similarity is quite natural since there is an intimate relation of supersymmetry to 
the Lorentz symmetries of space-time shown, for instance, by the fact that the anticommutators 
of the supersymmetry charges give translations. 

In the framework of quantum held theory the symmetries can be exploited by the corre¬ 
sponding Ward-Takahashi (WT) identities. Well studied examples are the global axial (chiral) 
symmetries in lattice QCD |59|. In case of SYM theory this way of realizing supersymmetry has 
been first considered by Curci and Veneziano |35| . At zero gaugino mass both supersymmetry 
and anomalous chiral symmetry has to be manifested in the corresponding WT identities. 

Before going to the lattice formulation let us describe the WT identities of supersymmetry 
in the continuum at a somewhat formal, non-rigorous level (for a more rigorous treatment see 
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(91) 


0). The corresponding supercurrent S tl has been introduced in (|5|): 

Sn(x) = -F p a r (x)a pr 7 M A a (x) = cr pT 7 M Tr (F^(ar)A(a:)) , 

where the held strength tensor matrix F )W {x) is dehned in (j7[) and the gaugino held matrix is 
A(;r) = A a (x)T a . Since the regularization breaks supersymmetry the simple conservation low 
d/jS^x) = 0 can only be recovered in the process of renormalization if the possible operator 
mixings are properly taken into account. The analysis of the operator mixings [J49] leads to the 
following conjecture about the form of SUSY-WT identities: 


+ Z r d ll T ll {x)0(y )) = (m 0 - m){D s {x)0(y)) . 


(92) 


Here 0(y ) is any gauge invariant local operator at point y and T^(x) is another dimension 
operator besides S^x), namely 


2 i 

T P (Y) = —7 v Tr (F^(x) A(x)) . 


(93) 


(94) 


The operator Dg(x ) appearing on the right hand side of (|92D is dehned as 

D s (x ) = ^cr pt/ Tr (F flu (x)X(x)) . 

Zg and Zt are multiplicative renormalization factors, m 0 is the bare gaugino mass and fh 
represents the additive renormalization of the bare mass. The form of ([)2j) is valid only for x j- y 
because “contact terms” have been omitted. The consequence of (^) is that the (renormalized) 
gaugino mass vanishes if m 0 — fh = 0 and the renormalized supercurrent can be dehned as 


S M (x) = ZgS^x) + ZtT^x) . 


(95) 


6.1 Lattice formulation 

The arguments leading to the SUSY-WT identity (|^) can be followed step by step in lattice 
regularization. First one has to define the supersymmetry transformations. Since SUSY is 
broken on the lattice the lattice action will not be invariant with respect to these transforma¬ 
tions. There is a freedom in the definition which is only restricted by the requirement that 
in the continuum limit the inhnitesimal transformations in (^') have to be reproduced. The 
supersymmetry breaking terms, which should vanish as O(a) in the continuum limit a —> 0, 
can be minimized by an appropriate choice of the irrelevant parts. An important requirement 
is to maintain as many discrete lattice symmetries as possible. In particular the parity V, time 
reversal T and charge conjugation C transformations should commute with the lattice SUSY 
transformations. 

A simple choice fulfilling these requirement has been introduced in [Q. For these defi¬ 
nitions let us change the conventions in the lattice action for gauginos (|43| ) by changing the 
normalization of the gaugino held according to 

A x —► y/m 0 + 4 A* (96) 
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and introduce the bare lattice mass tuq by 


K = 


(97) 


2 (mo + 4) 

Using the gaugino held matrix A^ = A“T a we have instead of (|43|) the fermionic action 

S, = E Tr { (mo + 4) XX, - i E [WM1 + 7e) } . (98) 

a: Z p= 1 

(Note that compared to ref. [Q U X4l is replaced here by U^.) This notation is close to the 
continuum conventions whereas (|43p is practical for numerical simulations. For the dehnition of 
the SUSY transformations we need an appropriately defined held strength tensor on the lattice 
which we denote by P^(x)\ 


P^w(x) = 


— [ u i^( x ) ~ U^(x) 


8i 9~i 

U\p v {x) = , 

U2[iv{x ) = U xfl U x —0-\-fi t uU x —0 t ijU x _p ll , 

&3/in(x) = , 

U^v^x) = Ux-ix,iiU x _^ y U x+ 0_^ fJI U xv . 


(99) 


This transforms under parity and time reversal in the same way as F )IV does in the continuum. 
Using P^v(x) one can define the inhnitesimal SUSY transformations on the lattice by 


5,U X „. = 


W 

J £ u xn — 2 L 

2 L 


X x + ^^A X+P.U, 


X[l 


s e ul = + 


S(Xx - 2 s " 1 "" 


pyPpv^x) . 


( 100 ) 


Here £ and £ are Grassmannian parameters satisfying a Majorana condition like (|2|). 

After defining the lattice SUSY transformations (|100|) one can derive the SUSY WT- 


identities by a standard procedure [^, |26|. The corresponding lattice supercurrent Sp(x) can 
be identified as 


1 


^fi{x) „ "y ] &Tr Pp a {x)U x ^\ x ^jJJ x ^ + Pp a {x + pyU X p\ x U x ^ 


( 101 ) 


pa 


The spinorial density multiplying the gaugino mass turns out to be 

Ds(x) = H<V Tr [TVMAJ ■ 


( 102 ) 


pa 
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The other supercurrent T^(x) which is mixed to SJx) can be defined as 


Tp{x) = 7^ Tr \P^v{x)Ul + P^(x + fi)U x ,\ x Ul 


Xfi 


(103) 


In the continuum limit, apart from differences in normalization conventions, these supercurrents 


go over to their continuum conterparts in (91)-( 

Using these definitions the SUSY WT-identity can be written as I0E3 0 : 

Z S (<S„(x)0(y)) + Z t (V‘T„(x) 0( 9 )) = (mo - m){D s (x)0(y)) + 0(a) . 


(104) 


Here the lattice (backward) derivative is defined as Vj)/(: x) = f{x) — f(x — ji). The gauge 
invariant function 0(y) at point y is assumed to be sufficiently far away from x in such a way 
that no common points with the expressions defined at x occur. In this way additional contact 
terms are avoided. 

The SUSY WT-identity as it stands in (|104| ) is only valid for gauge invariant functions 0(y) 
of lattice fields. In case of gauge non-invariance several additional terms appear. Examples of 
such anomalous terms emerge, for instance, in lattice perturbation theory calculations [[54], |55| 
The reason behind these complications is the conflict of gauge fixing with supersymmetry. 


As already remarked above, the form of the lattice currents in ( |101| )- (|103| ) is not unique 
because terms of order 0(a ) can be added. Locality and simplicity of the expressions is always 
welcome. This leads to the alternative definitions of the “local” supercurrents 


S l( x ) = -H<V 7 ax Tr [■ p P ai x )K} 


pa 


and 


hM=2y 7 „ Trl/VMAJ 


In contrast to these we can call the supercurrents in ( |101| ) and ( |103| ) as “point-split”: 

s»(x) = S„(x) , lyw = T„(x) . 


(105) 


(106) 


(107) 


In addition to the ambiguity of supercurrents the lattice approximation of the derivatives 
can also be varied. For instance, for the local supercurrents in (|105|) and (|106|) the backward 


difference Vqu may be replaced by a “symmetric difference” V®/(x) = |(/(x + ji) — f(x — fi)). 
All these variations are irrelevant in the continuum limit but for non-zero lattice spacing a/0 
the O(a) lattice artefacts in (|104|) can be rather different. 


6.2 Numerical results 

The DESY-Miinster-Roma Collaboration studied the SUSY WT-identities in numerical simu¬ 


lations H|, ffij. Omitting 0(a ) terms and dividing by Zs the SUSY WT-identity in (|104|) can 
be brought to the form ]27| 


(VjA^lOfe)) + ^(V%(x)0(y)) = (lf ‘° m) (Os(^)0(y)) . 

Zjq 


(108) 
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Figure 9: Results for the supercurrent renormalization factors from 
On the horizontal axis the inverse hopping parameter values 
are shown (l/n= 1/K). 


This has to be valid for every gauge invariant function O(y) which is defined in a point y in 
such a way that there are no common points with the functions defined at x. Considering (|108|) 
with different 0(y ) a system of linear equations is obtained for the two unknowns 


Z TS = 


^5 


_ {mo 

ZMS = - 7 


m 


(109) 


(Note that Zms is proportional to the renormalized gaugino mass uir = (mo — rh)/Z m .) Since 
there are very many different possible C?(i/)’s, the expectation that ( |108|) has, in the continuum 
limit, a unique solution pair ( Zts , z ms ) is highly non-trivial. Its numerical investigation can 
strongly support (or perhaps contradict) our expectations about the realization of supersym¬ 
metry in SYM theory. 

For obtaining non-zero expectation values the gauge invariant functions O(y) appearing in 
( |109| ) have to be chosen appropriately. Roughly speaking they should have quantum numbers 
similar to the spinorial density Ds(x). The intermediate states with important contributions 
are composite states of a gluino (= gaugino) and a gluon as, for instance, the gluino-glueball 
X investigated in section |0|. In order to obtain a good signal to noise ratio one can sum 
over three ^-coordinates (xi,X 2 ,x 3 ) keeping the fourth time-coordinate fixed. This projects 
out intermediate states with zero spatial momentum. One can also apply optimized smearing 
techniques in the timeslice containing y similarly to correlators for obtaining the masses of 
gluino-glueballs. 

The numerical simulations of the DESY-Miinster-Roma Collaboration j[ff]| has been per¬ 
formed on a 12 3 • 24 lattice at /3 = 2.3 for three different values of the bare gaugino mass given 
by K = 0.1925, 0.1940, 0.1950. These parameters are well suited for a first study because they 
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1/k 1/k 

Figure 10: Results for the renormalized gaugino mass from m- On 

the horizontal axis the invers hopping parameter values are shown 
(1/k = 1/K). The dashed lines are linear extrapolations for van¬ 
ishing gaugino mass. 


are reasonably close to the continuum limit and since the WT-identities hold in any volume 
no finite volume effects are disturbing. The values of the scale parameter r 0 |50], which give 
information about the physical volume size, are shown in table |j. A summary of the resulting 
fits to zts and zms is contained in figures |9] and |l0j. (For more details see 

The renormalization factors Zs and Zt depend on the definition of the lattice supercurrents. 
As figure |9] shows, the admixture of T ;i (x) is negligible for the point-split supercurrent S/ s (x) 
but non-zero for the local supercurrent S/f x). zts is within errors independent of the bare mass 
given by K. zms is proportional to the renormalized mass and tends to zero if K approaches 
the critical value K cr . 

These results are consistent with expectations, in particular with the realization of super- 
symmetry in a Yang-Mills theory with massless Majorana fermion in the adjoint representation. 
For a stronger conclusion the continuum limit /3 —> oo has to be investigated in future simula¬ 
tions. 


Table 2: The values of r 0 /a for the sets of configurations considered in 
ratio L/r o is also shown, where L is the spatial lattice size. 


The 


K 

r o/ a 

L/r 0 

0.1925 

6.71(19) 

1.79(5) 

0.194 

7.37(30) 

1.63(7) 

0.1955 

7.98(48) 

1.50(9) 
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7 Outlook 

The numerical Monte Carlo simulation of supersymmetric Yang-Mills theory is feasible with 
presently available computational resources and with known algorithms. The total computa¬ 
tional effort of the DESY-Munster-Roma Collaboration has been of the order of 10 Gflops • 
year ~ 3 • 10 1 ' floating point operations. This is much less than the amount of computer time 
devoted up to now to QCD. In fact, to simulate SYM theory is substantially easier than QCD 
because the “flavour number” of fermions Nf = | is smaller. In addition, in contrast to QCD, 
in SYM there are no almost massless pseudoscalar mesons expected. The easiness and the 
considerable theoretical interest make the Monte Carlo simulation of SYM theory a promising 
future subject. 

The next steps in improving the present simulations are obvious: 

• clarification of the nature of the transition at vanishing gaugino mass by studying its 
behaviour on larger lattice volumes; 

• investigation of the spectrum of bound states in supermultiplets in sufficiently large vol¬ 
umes and closer to the continuum limit; 

• study of the continuum limit of the supersymmetric Ward-Takahashi identities. 

Once these basic questions are sufficiently clarified one can certainly find other more detailed 
questions which will help understanding many aspects of the dynamics of supersymmetric gauge 
field theories. 
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